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DIFFERENTIAL  FLATNESS  AND  COOPERATIVE  TRACKING  IN  THE 

LORENZ  SYSTEM* 

LUIS  G.  CRESPO* 

Abstract.  In  this  paper  the  control  of  the  Lorenz  system  for  both  stabilization  and  tracking  problems 
is  studied  via  feedback  linearization  and  differential  flatness.  By  using  the  Rayleigh  number  as  the  control 
-only  variable  physically  tunable-  a  barrier  in  the  controllability  of  the  system  is  incidentally  imposed.  This  is 
reflected  in  the  appearance  of  a  singularity  in  the  state  transformation.  Composite  controllers  that  overcome 
this  difficulty  are  designed  and  evaluated.  The  transition  through  the  manifold  defined  by  such  a  singularity 
is  achieved  by  inducing  a  chaotic  response  within  a  boundary  layer  that  contains  it.  Outside  this  region,  a 
conventional  feedback  nonlinear  control  is  applied.  In  this  fashion,  the  authority  of  the  control  is  enlarged 
to  the  whole  state  space  and  the  need  for  high  control  efforts  is  mitigated.  In  addition,  the  differential 
parametrization  of  the  problem  is  used  to  track  nonlinear  functions  of  one  state  variable  (single  tracking) 
as  well  as  several  state  variables  (cooperative  tracking).  Control  tasks  that  lead  to  integrable  and  non- 
integrable  differential  equations  for  the  nominal  flat  output  in  steady-state  are  considered.  In  particular,  a 
novel  numerical  strategy  to  deal  with  the  non-integrable  case  is  proposed.  Numerical  results  validate  very 
well  the  control  design. 

Key  words,  feedback  linearization,  differential  flatness,  Lorenz  system,  cooperative  tracking,  non-linear 
control 

Subject  classification.  Applied  and  Numerical  Mathematics 

1.  Introduction.  The  analysis  and  control  of  chaotic  systems  have  attracted  considerable  attention 
in  recent  years.  A  chaotic  system  is  strongly  sensitive  to  small  changes  in  the  initial  conditions.  Such  a 
behavior  can  be  beneficial  or  detrimental  depending  upon  the  system  and  the  objective  under  investigation. 
In  this  paper  we  will  make  use  of  the  chaotic  response  of  the  Lorenz  system  to  enhance  the  performance  and 
controllability  of  feedback  linearization  based  controllers. 

The  Lorenz  system  is  a  simplified  model  of  a  thermally  driven  fluid  convection  system  between  parallel 
plates.  Depending  on  the  system  parameters,  such  a  system  exhibits  a  rich  spectrum  of  responses.  The 
control  of  the  Lorenz  system  has  been  studied  by  several  researchers  in  recent  years.  Vincent  and  Yu  [10] 
proposed  a  bang-bang  optimal  control  for  stabilizing  the  unstable  equilibrium  points  of  the  system.  Gao  and 
colleagues  [5]  studied  the  nonlinear  feedback  control  based  on  state  space  design.  Chen  and  Liu  [2],  Talwar 
and  Namachchivaya  [9]  and  Alvarez-Gallegos  [1]  investigated  the  nonlinear  regulation  of  the  Lorenz  system 
by  using  the  feedback  linearization  techniques  with  different  control  structures  and  control  objectives. 

From  the  physics  point  of  view,  it  is  natural  to  select  the  Rayleigh  number  as  the  control  variable. 
However,  such  a  strategy  makes  the  system  uncontrollable  at  the  plane  x  =  0.  This  control  variable  was 
used  by  the  author  in  [3].  When  feedback  linearization  techniques  are  applied,  this  feature  appears  as  a 
singularity  in  the  state  transformation.  When  the  feedback  linearization  techniques  are  applied,  this  feature 
also  appears  as  a  singularity  in  the  state  transformation.  The  controllability  of  the  system  is  limited  to  half 
of  the  state  space.  The  singularity  is  also  responsible  for  extremely  high  control  efforts  in  its  vicinity.  In  this 
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paper,  we  propose  composite  controllers  to  overcome  this  limitation.  Within  a  boundary  layer  that  contains 
such  a  singularity,  the  chaotic  response  of  the  system  and  its  response  to  step  inputs  are  used  to  drive  the 
system  through.  Furthermore,  the  differential  flatness  of  the  system  is  used  to  aim  for  control  objectives 
that  do  and  do  not  admit  a  closed  form  expression  for  the  corresponding  flat  output  of  the  system. 

This  paper  is  organized  as  follows.  A  brief  description  and  analysis  of  the  Lorenz  system  are  presented 
in  Section  2.  Section  3  discusses  the  feedback  linearization  and  its  application  to  the  Lorenz  system.  In 
Section  4,  composite  controllers  are  designed  and  evaluated  for  both  stabilization  and  tracking  problems. 

2.  The  Lorenz  System.  The  Lorenz  model  is  obtained  from  studying  a  fluid  layer  heated  from  below 
and  cooled  from  above  such  that  a  temperature  difference  is  established  across  it.  The  convection  motion 
is  described  by  the  Navier-Stokes  equations.  Taking  Fourier  expansion  of  these  equations  along  two  spa¬ 
tial  directions  and  truncating  the  remaining  expressions  to  retain  only  three  modes  leads  to  the  following 
simplified  model 


x  =  cr(y  —  x), 

y  =  px-y-xz ,  (2.1) 

i  =  -j3z  +  xy , 


where  a,  p  and  3  are  real  parameters  denoting  the  Prandtl  number,  the  Rayleigh  number  and  a  geometric 
factor,  respectively.  The  state  variables  x,  y  and  z  represent  measures  of  fluid  velocities  and  the  spatial 
temperature  distribution  in  the  fluid  layer  under  gravity.  From  the  physical  point  of  view,  the  Rayleigh 
number  p  can  be  easily  manipulated  by  changing  the  heat  transfer  to  the  fluid  from  below.  This  parameter 
will  be  treated  as  the  control  variable.  We  denote  u  =  p. 

For  the  brevity  of  discussion,  we  assume  that  the  control  converges  to  a  constant  value  in  steady  state 
denoted  as  uss  =  lim^oott.  The  singular  points  of  the  system  can  be  parametrized  with  uss .  The  locus  of 
these  points  is  given  by  a  family  of  curves 

XI  =  [0,0,0]T,  xy  =  [±V/3(Uss  -  1), ±V/?(u„  - !),«..  - 1  ]T-  (2.2) 

Linearization  about  X{  leads  to  characteristic  equation  A3  +  AA3  +  BX  +  C  -  0,  where  A  =  1  +  ft  +  <t, 
B  =  j3(a  + 1)  +  a  —  uss  and  C  =  /3(<j  —  uss )•  Linearization  about  X leads  to  A  =  1  +  /?  +  cr,  B  =  0(uss  +  a) 
and  C  —  2 (3a(uss  -  1). 

The  stability  analysis  leads  to  the  following  observations.  When  uss  <  1  the  origin  of  the  system  is  a 
stable  equilibrium  point.  When  1  <  uss  <  P  =  cr(a  +  p  +  3)/(<r  —  0  —  1),  X±  is  unstable  and  X^3  are  stable. 
When  uss  >  p  there  are  no  stable  equilibrium  points  and  the  system  reaches  a  chaotic  regime. 

At  u8S  =  1  a  pitchfork  bifurcation  takes  place  at  X*,  while  for  uss  —  p  a  subcritical  Hopf  bifurcation 
occurs  at  X|?3.  For  uss  >  p,  the  system  is  driven  by  repulsions  exclusively  while  the  trajectories  are  confined 
to  a  region  of  finite  volume  forming  a  strange  attractor.  The  response  on  the  attractor  is  chaotic.  For 
additional  information,  the  reader  can  refer  to  [8].  In  the  numerical  simulations,  the  parameters  to  be  used 
are  cr  =  10  and  j3  =  8/3. 

3.  Feedback  State  Linearization. 

3.1.  Background.  Consider  the  single  input  system 

x  =  f(x)  +  g(x)it,  (3.1) 
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where  x  G  Rn  is  the  state  vector,  u€  R  is  the  control  and  f ,  g  :  Rn-»  Rn  are  sufficiently  smooth  nonlinear 
functions  of  their  arguments.  The  Lie  derivative  of  g(x)  with  respect  to  the  vector  field  f(x)  is  defined  as 

adfg(x)  =  [f,  ad£_1g](x)  for  k  >  1, 

a<$g(x)  =  g(x),  (3.2) 

where  [X,  Y]  =  VY  *  X-VX  •  Y  is  the  Lie  bracket  of  the  vector  fields  X  and  Y.  According  to  [6],  there 
exists  a  real  value  function  A(x),  defined  in  a  neighborhood  U(xq)  of  xo  such  that 


Le  A(x)  =  Ladtg  A(x)  =  ...  =  Ladr  2gA(x)  =  0,  (3.3) 

^-‘^(xo)  /  0, 

where  LgA(x)  =  VA(x)  •  g  denotes  the  Lie  derivative  of  the  real- value  function  A(x)  with  respect  to  the 
vector  field  g,  if  (i)  the  matrix 

C  =  [g(x0),  adfg(xo),  adfg(xo), ...,  adf_1g(x0)],  (3.4) 

has  rank  n  and  (ii)  the  distribution  D  =span{g(x),adfg(x),adfg(x),  ...,adf  _1g(x)}  is  involutive.  Further- 
more,  there  exists  a  transformation  z  =  <f>(x)  in  U(xq)  such  that 

$(x)  =  [A(x),  If  A(x), ...,  If  _1A(x)]t,  (3.5) 

dn_1A(x).T 


=  [A(x),A(x),...,- 


dtn~l 


The  system  (3.1)  is  transformed  into  a  chain  of  integrators 

ii  =  z2,  z2  =  z3,  ■■■  ,  i„_i  =  zn,  zn  =  o(x)  +  b(x)u  =  v, 


(3.6) 


v(x)  =  a(x)  +  b(x)u  —  IgIf  1A(x)  +  If  A(x)u,  (3.7) 

where  Zj  =  If-1  A(x)  (*  =  1, . . . ,  n).  The  control  for  this  linear  system  can  be  designed  by  full  state  feedback 
and  pole  placement  techniques.  For  example,  the  control  v  can  be  taken  as: 

u($"1(z))--Kz,  (3.8) 

where  the  feedback  gains  K  are  chosen  to  place  the  closed  loop  poles  at  the  desired  locations  in  the  left  hand 
side  of  the  complex  plane.  For  tracking  problems,  the  feedback  signal  will  be  the  tracking  error  z^  —  z  where 
z d  is  a  pre-specified  reference  signal  in  the  z- domain. 

The  control  in  the  physical  domain  is  then  given  by: 

u(x)  =  [n($_1(z))  -  o(x)]  /b(x).  (3.9) 

It  should  be  noted  that  u(x)  becomes  unbounded  when  b(x)  -»  0. 

3.2.  Feedback  Linearization  of  the  Lorenz  System.  For  the  Lorenz  system,  f  =  [a(y  —  x),  -y  — 
xZj  —fiz  4-  xy]T  and  g  =  [0,  x,  0]T.  After  some  manipulations  we  obtain 

adg(adfg)  =  [0, 2a x,  0]T, 

adfg  =  [— ax,a(y  —  x)  +  x,  —  x2]T :  (3.10) 

a(a  -  l)x  -  2a2y 

- x 3  -  2axy  +  (a  -  l)2x  +  a(l  -  a)y 
(3<j  -  1  —  p)x2  —  2axy 
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By  evaluating  the  rank  of  the  C  matrix  in  Equation  (3.4),  we  find  that  the  rank  is  3  except  when  the 
transformation  is  singular  at  x  —  0  or  a  =  (3/2.  At  the  singularity  at  x  =  0,  the  system  (2.1)  is  completely 
insensitive  to  the  control.  This  restriction  imposes  an  unavoidable  barrier  in  the  controllability  of  the  system, 
leaving  just  half  of  the  state  space  at  the  disposal  of  the  control.  Which  half  the  system  stays  in  depends  on 
the  initial  condition.  The  closer  the  system  gets  to  this  singular  plane,  the  higher  the  control  effort  will  be, 
approaching  positive  or  negative  infinity  as  x  approaches  to  zero. 

From  Equation  (3.3),  the  function  A(x)  ,  known  as  the  flat  output  of  the  system  [4,  7],  satisfies 


dX 


DX  n 

3  77-  =  0, 

Oy 

dX  2dX 

_  =  0. 


(3.11) 


Solving  this  set  of  equations,  we  obtain  A  —  x2  /  2  -  az  4-  fc,  where  k  is  the  integration  constant.  Equation 
(3.5)  leads  to  the  state  transformation 


z=[zi,z2,z3}t  =  $(x)  (3.12) 

=[x2 /2  —  az  4-  fc,  a(fiz  —  x 2),  ayxy  4*  2<j2x2  —  afi2z]T , 


where  7  =  fi  —  2a.  The  inverse  transformation  is  given  by 

x=[x,y,z]T  =  $-1(z)  (3.13) 

_  +  /2(/fe  +  +  2a/3zi  +  (/?  +  2a )z2  +  Z3  -  2 a0k  2ae  +  z2 

V  7  ’  a^/2'-f(^£  + z2) 

where  e  =  z\  -  k.  The  transformed  dynamic  system  takes  the  form  of  Equation  (3.6)  with  n  =  3  and  a(x) 
and  6(x)  in  Equation  (3.7)  given  by 


a(x)  =  -ayx2z  -  4a3x2  -  a( 7  4-  ft2  +  afi  -  6a2)xy  4-  a2yy2  +  afi3z,  (3.14) 

b(x)  =  ayx2. 


By  using  the  solution  of  the  PDE  (3.11)  and  the  system  equation  (2.1),  we  have 

A  =  x2/2  -  <7 z  +  k,  (3.15) 

A  =  a(fiz  -  x 2), 

A  =  ayxy  +  2a2x2  —  fi2az. 


The  state  variables  and  the  control  can  be  differentially  parametrized  using  the  flat  output. 
x  =  ±y/2x] 7, 

y  =  ±(A  4-  03  +  2cr)A  4-  2afiX  -  2a/e)/(av/2^),  (3.16) 

z  =  (fiX  4-  2a fiX  -  2ak)/(afiy ), 

u  ~  - — \ — ^{~2fiyxv  -  2fiy(fix  4-  x  +  2ca  -  crA;)A  -  2fi(4kc  -  2 c2  -  2 k2  4-  2ayx  + 

4afiyxz 

fiyx  4-  2afiyc  —  yafik)X  —  4afi(fiyx  +  4kc  —  2c2  —  2ft2)A  —  Sake2  +  Aafiykx  —  8  ak3  + 

16afe2c  4~  fi2y(fi  +  2cr)A2  4-  2fiy(fi  4~  (j)AA  H-  fiyX2  4"  2afi2yXX  4~  2afi3yXX}, 

where  c  =  A  4-  /3A,  x  =  c  -  k  and  v  =  X .  Inequality  range  constrains  on  u  can  be  imposed  by  designing 
composite  controllers.  This  is  explained  in  the  next  section. 
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4.  Composite  Control.  From  the  above  discussions,  we  know  that  a  feedback  linearization  based 
control  is  not  able  to  drive  the  system  across  the  singularity  imposed  by  the  transformation.  This  fact 
restricts  the  controllability  of  the  system  to  a  portion  of  the  state  space,  half  in  this  case,  and  leads  to 
extremely  high  control  effort  in  the  vicinity  of  the  singularity.  In  this  section  we  propose  a  composite 
controller  to  overcome  such  a  difficulty. 

Denote  the  hyper-plane  that  makes  the  state  transformation  singular  as  g(x)  =0.  Define  the  composite 
controller  given  by: 

n(x)  =  {Ml  if|x-p(x)|>*  (4.1) 

I  U2  otherwise 

where  S  defines  the  thickness  of  a  boundary  layer  about  the  singular  plane,  U\  is  the  control  expression 
resulting  from  using  a  conventional  control  method  and  U2  is  the  amplitude  of  a  step  input  set  according  to 
the  particular  control  objectives. 

The  condition  for  applying  u\  can  be  replaced  by  |ui|  <  U ,  provided  that  control  saturation  occurs  in 
the  vicinity  of  the  boundary  layer.  In  this  paper,  g(x)  =  0  and  U2  is  a  step  input  that  induces  a  chaotic 
response  within  the  boundary  layer  i.e.  U2  >  p. 

The  non-empty  intersection  of  the  attractor  and  the  two  hyper-planes  defined  by  |x— </(x)|  =  S  guarantees 
that  the  crossing  of  the  boundary  layer  occurs.  This  can  be  proved  as  follows.  Let’s  call  hi(x)  for  i  —  1,2 
these  two  planes.  Assume  that  (i)  the  chaotic  response  is  moving  within  a  strange  attractor  whose  state 
space  location  is  given  by  A\  and  that  (ii)  Ai  ft  h*(x)  ^  0.  The  crossing  of  the  hyper-plane  ft*(x)  will  not 
occur  iff  there  exist  a  subset  A2  C  A\  such  that  A2  fl  h{(x)  =  0  for  all  times  i.e.  A2  is  a  strange  attractor 
by  itself.  This  implication  violates  the  irreducibility  property  of  A\  then  A2  can  not  exist.  The  need  for 
Ai  n  hi(x)  ^  0  imposes  bounds  to  6  from  above. 

Notice  that  in  this  scheme  the  control  does  not  have  the  authority  to  manipulate  the  transient  part  of 
the  transition  from  one  side  of  the  boundary  layer  to  the  other  one.  For  some  states  several  crossings  of  g(x) 
might  occur  before  the  system  leaves  the  boundary  layer.  Such  behavior  is  clearly  undesirable.  Refinements 
and  improvements  of  the  control  within  the  boundary  layer  can  be  achieved  by  taking  into  consideration 
the  control  objective  and  the  state  of  the  system.  This  practice  however  was  not  implemented  in  this  study. 
Once  the  system  leaves  the  boundary  layer  u\  is  applied  and  the  stabilization/tracking  is  achieved.  In  this 
paper  ui  is  given  by  u  in  Equation  (3.16).  Due  to  the  structure  of  the  controller,  global  uniform  asymptotic 
stability  about  the  reference  signal  r(x(t))  is  achieved  at  the  locations  where  the  intersection  of  r(x(t))  and 
the  boundary  layer  is  an  empty  set. 

For  stabilization,  the  linear  system  given  by  Equations  (3.6),  (3.7)  and  (3.14)  can  be  controlled  by  pole 
placement  techniques.  Taking  the  feedback  control  v  as: 

v  -  aiZi  +  a2z2  +  0:3^3,  (4*2) 

where  the  feedback  gains  ai  are  chosen  to  place  the  closed  loop  poles  in  the  left  hand  side  of  the  complex  plane. 
On  the  original  state  variables,  the  control  can  be  obtained  after  substituting  Equations  (3.12)  and  (3.8)  into 
Equation  (3.7).  Since  this  procedure  stabilizes  z,  the  steady-state  values  in  the  x-domain  can  be  controlled  by 
manipulating  k  according  to  Equation  (3.13).  The  reader  must  notice  that  stabilization  about  the  origin  using 
feedback  linearization  requires  infinite  control  effort,  i.e.  limx_»o  u(x)  =  lima._+o(v(x)  -  a(x))/b(x)  =  ±oo. 
However,  from  the  stability  analysis  we  know  that  any  control  satisfying  0  <  uss  <  1  will  drive  the  system 
to  the  origin  from  a  given  initial  condition. 
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4.1.  Single  State  Tracking  Control.  In  this  section  we  present  examples  for  stabilization  and  track¬ 
ing  control  problems  that  involve  a  single  state.  The  tracking  signal  Xd(t)  =  a  +  6sin(£)  for  the  state  variable 
x(t)  is  considered  here.  The  flat  output  corresponding  to  the  tracking  signal  Xd(t)  is  given  by 

A d(t)  =  (7/2/?) (62/ 2  +  a2)  +  (7 ab/(l32  +  1))  {/ 3 sin(t)  -  cos(t)}  (4.3) 

~  (b2'y/(2/32  4-  8)){(j8/2)  cos(2t)  -  sin(2i)} 

The  difference  between  the  system  flat  output,  Ad  -  A  =  Zd  -  z\  is  the  tracking  error.  We  take  the  full  state 
feedback  control  v  for  the  tracking  problem  as  follows 

v  =  Ad  +  /fi (Ad  -  A)  + 1^2 (Ad  -  A)  +  I<3(Xd  -  A)  (4.4) 

The  control  gains  K{  are  selected  such  that  the  tracking  error  vanishes  exponentially.  It  should  be  noted 
that  the  tracking  control  is  designed  in  the  transformed  space  z,  and  therefore  is  indirect  for  x. 

The  state  variables  and  the  nominal  control  obtained  from  Equation  (3.16)  are  shown  in  Figure  5.1  for 
a  —  5  and  6  =  6.  The  control  grows  unbounded  near  the  singularity.  Time  evolutions  for  different  Xd(t) 
are  shown  in  Figures  5.2  and  5.3  using  5  =  0.1  and  S  =  0.2  respectively.  In  all  the  cases  the  control  is 
activated  after  30  seconds  of  chaotic  regime.  In  the  first  case,  the  system  does  not  reach  the  boundary  layer 
and  perfect  tracking  is  achieved  after  a  short  transient.  If  the  control  is  activated  when  the  singular  plane 
is  between  the  state  of  the  system  and  xd(t),  the  system  would  reach  and  cross  the  boundary  layer  before 
settling  down.  This  control  was  designed  such  that  u2  is  applied  when  u\  <  6.  In  this  fashion  the  control 
range  constraint  u  >  0  is  imposed.  In  the  second  case,  the  desired  trajectory  crosses  repeatedly  the  singular 
plane  forcing  the  system  to  reach  the  boundary  layer  several  times.  The  effect  of  not  applying  the  nominal 
control  is  slightly  noticeable.  The  reader  must  notice  that  both  the  x  and  y  states  behave  similarly  once  the 
system  reaches  the  strange  attractor. 

Figure  5.4  shows  the  results  of  the  composite  controller  with  a  wider  boundary  layer  5  =  0.5.  Recall 
that  the  control  does  not  have  authority  within  the  boundary  layer  and  relies  on  the  chaotic  behavior  of 
the  system  to  cross  it.  As  expected,  the  increase  of  6  has  a  detrimental  effect  on  the  control  performance. 
For  even  higher  values  of  (5,  the  system  might  be  trapped  in  the  boundary  layer  while  the  desired  trajectory 
completes  half  a  cycle. 

4.2.  Cooperative  Tracking  Control.  Now  we  use  the  differential  parametrization  of  the  states  and 
control  to  aim  for  tracking  objectives  that  involve  combinations  of  the  states.  The  problem  statement  is  as 
follows.  Find  ui  in  Equation  (4.1)  such  that  the  system  is  driven  to  the  manifold  defined  by  h(x,u,t)  —  0 
from  any  initial  condition.  In  this  problem,  the  system  would  track  signals  that  imply  cooperative  relations 
among  the  states,  being  the  tracking  of  a  trajectory  of  a  particular  state  a  particular  case.  Again,  the 
composite  structure  of  the  control  enable  to  achieve  global  stability  about  the  desired  tracking  function.  In 
the  examples  to  come  we  will  take  u2  —  30  >  p  within  the  boundary  layer. 

Depending  upon  the  tracking  objective,  the  equation  for  the  desired  flat  output  can  be  integrable  or 
non-integrable.  Both  cases  are  considered  next. 

4.2.1.  Problems  with  closed  form  solution  for  A<*.  Once  the  control  objective  /i(x,u,£)  —  0  is  set, 
the  Equations  (3.16)  lead  to  an  ODE  for  the  desired  flat  output  A d.  Such  an  equation  is  in  general  non-linear. 
In  this  section  we  study  problems  where  a  closed  form  expression  for  the  steady  state  solution  can  be  found. 

As  the  first  example,  we  take  h(x,u,t)  —  x2  —  z  —  a  —  bsin(ut).  The  corresponding  differential  equation 
for  the  flat  output  is  given  by 

A  d  +  (2a^/(j))\d  -  007/ <{>  +  (baj/cj))  sin  (art), 


6 


where  (j)  ~  2a  —  1  and  £  =  f3  —  1 .  Solving  we  find 

A d  =  a7/(2£)  +  ??exp{  — 2£<t£/</>}  +  {25<t27£  sm(ut)  —  cobajficos^t)} /(4<t2£2  +  u?2<^2).  (4.5) 

For  our  purposes,  providing  that  the  transient  response  vanishes,  only  the  steady  state  component  is  needed. 
This  expression  along  with  Equations  (3.16)  and  (4.4)  fully  determine  the  control.  After  some  manipulations 
we  find  that  the  system  reaches  stationarity  at 

xsa  =  .  (4-6) 

^ ss  —  2  +  (a  — *  £)/£, 

Numerical  results  for  a  =  15  and  5  —  5  are  shown  in  Figure  5.5.  For  this  particular  case,  the  system  crosses 
the  boundary  layer  twice.  Figure  5.6  shows  the  numerical  results  for  a  =  10  and  5  =  0. 

As  a  second  case,  we  use  h(x)  =  x2  +  y2  4-  z2  —  R2.  This  problem  can  be  interpreted  as  the  stabilization 
of  the  system  about  the  surface  of  a  sphere.  The  corresponding  differential  equation  for  A<*,  not  shown  here, 
is  non-linear.  However,  because  the  homogeneous  solution  vanishes  with  time  we  can  solve  for  particular 
solution  and  use  it  in  Equation  (4.4)  to  calculate  the  control.  The  corresponding  flat  output  and  steady 
state  values  are: 


A d  =  7(±0-/3)/2, 

xss  =  [±v//J  {6  -  J),  ±v//3  (0  -  /?),  0  -  /?]T , 

v>ss  =  ±(l  +  0-/3), 

where  9  =  \J f32  +  R2.  Stationary  values  of  the  states  can  also  be  found  by  solving  for  the  values  of  k  in 
Equation  (3.13)  that  satisfy  h(x)  =  0  when  z  — >  0  (due  to  Equation  (3.8).  Numerical  results  for  R  =  2  are 
shown  in  Figure  (5.7).  In  the  case  shown,  the  boundary  layer  is  not  reached. 

4.2.2.  Problems  without  closed  form  solution  for  A^.  In  the  previous  section,  the  tracking  objec¬ 
tives  led  to  stable  ODEs  for  A d  whose  particular  solution  could  be  found  in  closed  form.  In  this  section,  we 
consider  problems  in  which  this  is  not  the  case. 

By  integrating  numerically  and  simultaneously  both  the  state  Equations  (2.1)  and  the  ordinary  differen¬ 
tial  equation  for  the  desired  flat  output  A<*,  tracking  can  be  achieved.  It  is  important  to  notice  that  realizable 
objectives  imply  stable  solutions  for  A <*.  The  tracking  error  dynamics,  set  by  Equation  (4.4),  makes  the  flat 
output  in  Equation  (3.15)  converge  to  the  steady  state  value  of  A^.  Tracking  is  achieved  once  the  transient 
for  both  the  real  and  the  desired  flat  outputs  vanish. 

Due  to  the  non-linear  character  of  the  equation  for  A <*,  the  steady  state  response  might  exhibit  dependence 
to  the  initial  conditions.  For  this  reason,  the  control  design  must  start  by  searching  for  the  initial  conditions 
in  A d  and  its  derivatives  that  lead  to  steady  state  trajectories  that  satisfy  the  desired  response  specifications. 
The  reader  must  notice  that  a  fixed  number  of  derivatives  of  A <*,  three  in  this  case,  are  needed  to  build  v. 
This  can  be  obtained  by  performing  additional  differentiations  on  the  ODE  for  Ad. 

As  an  example,  we  use  the  energy-like  expression  h(x,t)  =  x2  +  y2  +  mz  -  E(t).The  corresponding 
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differential  equation  for  Ad  is  given  by: 


0  =  -2EP^2a2(Xd  +  A d0  -  k)  +  2mj<rl3Xd  -  imkja2Xd  +  im-ya2 (3XdXd  + 

4mya2  ft2  Xd  +  4mk2ija2  +  2myaP2  XdXd  —  2m,yaPkXd  —  8m^a2PkXd  + 

24ft2  k2  o2  Xd  -  24a2  PkX2d  +  24k2a2f3Xd  -  48a2kp2XdXd  +  8a2pX3d  -  (4.7) 

8k3 a2  (i  +  4p3ya2X2d  +  -f/3X2d  +  2'y01XdXd  +  4yapXdXd  +  4'yap3XdXd  - 
8kja2pXd  +  4'fafi2X2d  +  4k2ja2p  -  4k^aPXd  -  4kyap2Xd  +  8ya2p2XdXd  + 

47<t2/3A^  -  8P2kja2Xd  +  4'y(T/32XdXd  +  -yp3X2d  +  8/34o2Xd  +  24p3a2X2dXd  - 
24  p3a2kX2d  +  24a2p2XdX2 

The  set  of  equations  to  be  integrated  numerically  is  given  by  Equations  (2.1)  and  (4.7),  where  u  is  given 
by  Equations  (4.1)  and  (4.4).  The  reader  must  notice  that  does  not  depend  on  the  states  and  that  both 
equations  are  coupled  via  u.  Numerical  simulations  with  E(t)  =  a  +  bsin(wt)  cos(wt)  are  shown  in  Figure 
(5.8). 


5.  Conclusions.  This  paper  studies  the  stabilization  and  tracking  control  of  the  Lorenz  system  using 
feedback  linearization  and  differential  flatness.  When  the  Rayleigh  number  is  used  as  the  control  variable, 
the  system  is  uncontrollable  in  a  manifold  of  the  state  space.  In  the  vicinity  of  such  a  singularity,  the  control 
demands  grow  unbounded.  Composite  controllers  that  use  feedback  linearization  and  the  system  response 
to  step  inputs  are  proposed  to  overcome  this  difficulty.  By  inducing  the  chaotic  response  within  a  boundary 
layer  which  contains  the  singular  plane,  the  transition  to  desired  states  is  achieved.  Such  controls  can  be 
used  not  only  to  enlarge  the  controllability  region  of  the  system  to  the  whole  state  space  but  also  to  mitigate 
high  control  demands.  Control  objectives  and  initial  conditions  that  imply  single  and  multiple  crossings  of 
the  boundary  layer  are  studied  in  the  examples.  In  addition,  tracking  control  problems  that  involve  single 
and  cooperative  relations  among  the  states  are  studied  using  the  differential  flatness  of  the  system.  Problems 
with  control  objectives  that  lead  to  integrable  and  non-integrable  differential  equations  for  the  desired  flat 
output  are  considered.  A  numerical  approach  in  which  the  state  equations  and  the  differential  equation  for 
the  nominal  flat  output  are  simultaneously  integrated  is  proposed  and  validated.  Numerical  simulations  led 
to  excellent  performances. 

Acknowledgments.  The  author  would  like  to  thank  Professor  Sunil  K.  Agrawal  from  University  of  Delaware 
for  his  feedback  and  enthusiasm. 


REFERENCES 

[1]  J.  Alvarez-Gallegos,  Nonlinear  regulation  of  a  Lorenz  system  by  feedback  linearization  techniques , 

Dynamics  and  Control,  4  (1994),  pp.  277-298. 

[2]  L.  Chen  and  Y.  Liu,  Control  of  the  Lorenz  chaos  by  the  exact  linearization ,  Applied  Mathematics  and 

Mechanics,  19  (1998),  pp.  67-73. 

[3]  L.  G.  CRESPO  and  J.  Q.  Sun,  Feedback  linearization  of  the  lorenz  system:  Stabilization  and  tracking , 

Submitted  to  the  Journal  of  Vibration  and  Control,  (2002). 

[4]  M.  Fliess,  J.  Levine,  P.  Martin,  and  P.  Rouchon,  Flatness  and  defect  of  non-linear  systems : 

Introductory  theory  and  examples ,  International  Journal  of  Control,  61  (1995),  pp.  1327-1361. 


8 


[5]  F.  Gao,  W.  Q.  Liu,  V.  Sreeram,  and  K.  L.  Teo,  Nonlinear  feedback  control  for  the  Lorenz  system , 

Dynamics  and  Control,  11  (2001),  pp.  57-69. 

[6]  A.  IsiDORl,  Nonlinear  Control  Systems ,  Springer- Verlag,  New  York,  1989. 

[7]  H.  Sira-Ramirez,  R.  Castro,  and  E.  Liceaga,  Liouvillian  systems  approach  for  the  trajectory 

planning-based  control  of  helicopter  models ,  International  Journal  of  Robust  and  Non-linear  Control, 
10  (2000),  pp.  301-320. 

[8]  C.  Sparrow,  The  Lorenz  Equations:  Bifurcations ,  Chaos  and  Strange  Attractors ,  Springer- Verlag,  New 

York,  1982. 

[9]  S.  Talwar  and  N.  S.  Namachchivaya,  Control  of  chaotic  systems:  Applications  to  the  Lorenz 

equations ,  Proceedings  of  Nonlinear  Vibrations,  ASME,  50  (1992). 

[10]  T.  L.  Vincent  and  J.  Yu,  Control  of  a  chaotic  system ,  Dynamics  and  Control,  1  (1991),  pp.  35-52. 


9 


Fig.  5.1.  Nominal  time  histories  of  u  (solid  line),  x  (dot-dashed  line),  y  (dotted  line),  and  z  (dashed  line)  when  the 
reference  trajectory  is  Xd(t)  =  5  +  6sin(u;£). 


Time 


Fig.  5.2.  Time  evolutions  of  the  state  variables  and  the  control  with  6  =  0.1,  a  =  8  and  b  =  5.  The  control  is  activated 
after  30s.  The  system  does  not  reach  the  boundary  layer . 
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Fig.  5.3.  Time  evolutions  of  the  state  variables  and  the  control  with  6  —  0.2,  a  =  0  and  b  =  10.  The  control  is  activated 
after  30s.  Sharp  peaks  in  the  control  curve  indicate  that  the  system  reaches  the  singularity. 
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Fig.  5.4.  Time  evolutions  of  the  state  variables  and  control  with  5  =  0.5,  a  =  0  and  b  =  10.  The  control  is  activated  after 
30s.  Sharp  peaks  in  the  control  curve  indicate  that  the  system  reaches  the  singularity. 
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Fig.  5.5.  Time  evolutions  of  the  state  variables  and  the  control  for  tracking  the  hyper-plane  /i(x,  t)  —  x2  —  z  —  15  —  5sin(a;t). 
The  following  conventions  are  used  xfdotted  line),  y (dot-dashed  line),  zfdashed  line),  ufsolid  line)  and  h(thick  line) 


Time 

Fig.  5.6.  Time  evolutions  of  the  state  variables  and  the  control  for  tracking  the  hyper-plane  h(x)  —  x2  —  z  —  10.  The 
following  conventions  are  used  xfdotted  line),  y (dot-dashed  line),  z (dashed  line),  u(solid  line)  and  hfthick  line). 


Fig.  5.7.  Time  evolutions  of  the  state  variables  and  the  control  for  tracking  the  hyper-plane  h(x)  =  ar  -f  2r  +  —  4.  The 

following  conventions  are  used  x  (dotted  line),  y (dot-dashed  line),  z  (dashed  line),  u (solid  line),  and  h  (thick  line) 


Fig.  5.8.  On  the  top,  time  evolutions  of  the  states  xfdotted  line),  yfdashed  line)  and  z(solid  line)  are  shown.  In  the 
middle,  time  evolutions  for  the  flat  output  and  its  derivatives  are  shown.  In  the  bottom,  time  evolutions  of  E(t)  for  both,  the 
desired  (solid  line)  and  the  real  (dotted  line)  flat  outputs  are  shown. 
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